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Abstract: Under the assumption of Minimal Flavour Violation (MFV), the Unitarity 
Triangle (UT) can be determined by using only angle measurements and tree-level observ- 
ables. In this respect, the most accurate quantities today available are sin 2(3^k s , arm 
\V u b\. Among the latter, \V u b\ is at present the quantity suffering from largest systematic 
uncertainties, given the discrepancy between the inclusive and the exclusive determinations. 
We show with a numerical fit how sensitively the MFV-UT determination depends on 
the choice of \V u b\. In addition, we focus on the implications of the inclusive value for 
\V u b\, which favors two non-SM like solutions in the 7j — rj plane. We study in detail the 
possibility of reproducing such solutions within the MFV MSSM. Our findings indicate 
that the case for the MFV MSSM is in this respect quite problematic, unless the non- 
perturbative parameters £ and Bk are significantly different from those obtained by lattice 
methods. As a byproduct, we point out that scenarios with 200 GeV < Ma % 500 GeV 
and tan /3 — 50 that predict a significant suppression for AM S in correlation with an 
enhancement for BR(i? s — ► /x + /i~) have to be fine-tuned in order not to violate the new 
combined bound on the latter decay mode from the CDF and D0 collaborations. Relatively 
large correlated effects can however still occur for negative values of fi and large values for 
Ma > 500 GeV, increasing with increasing tan f3 > 30. 
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1. Introduction 



Most of the calculable models of low-energy New Physics (NP) available at present tend 
to be invasive in the sector of flavour violation and to destroy the very specific pattern of 
flavour changing neutral current (FCNC) effects predicted by the SM. In the example of 
SUSY, this happens because of its soft sector. In fact, in absence of compelling symmetries 
restricting the form of the soft terms, the latter must be parameterized most generally. 
However, when calculating SUSY contributions to flavour observables with such generic 
soft terms, it is hard to hide simultaneously all their effects behind the already small SM 
predictions, which agree rather well with the experimental data. The soft terms parameter 
space that turns out to survive FCNC constraints looks then very fine tuned and this raises 
the SUSY flavour problem. 

A different, phenomenologically coherent approach, is to assume - on the basis of the 
success of the SM description of FCNC processes - that the flavour sector of any NP model 
maintains a natural mechanism of "near- flavour-conservation" Q. The latter would then 
allow small effects in exact analogy with what happens for the SM. 

This idea has been first elucidated as a meaningful requirement for models of NP 
at the EW scale by the authors of Refs. j|, [jj. It has then been formalized in || as 
a consistent effective field theory framework, called Minimal Flavour Violation (MFV). 
In this framework, the SM mechanism of near-flavour-conservation is extended to NP by 
means of the SM Yukawa couplings. Since the latter are in the SM the only sources of 
flavour breaking, in MFV they play the role of fundamental "building blocks" of flavour 
violation and every new source of flavour violation, entering the NP Lagrangian, becomes 
function of them. 

This approach is more general than the so-called constrained MFV (CMFV) §, |, § , 
in which one also imposes the same operator structure as in the SM. 

Recently, the MFV framework of || has been applied to the MSSM at low tan f3, in a 
detailed numerical study of meson mixings 0. 1 Thereby, the phenomenological differences 
between MFV and CMFV have been spelled out. It has also been shown how the imple- 
mentation of the MFV limit makes SUSY contributions to meson mixings naturally small, 
even for a SUSY scale of a few hundreds GeV. 

If NP is of MFV nature - as the lesson of 5-factories seems to hint - the findings of 
Ref. indicate that the search for NP effects in flavour observables is more challenging, 
but not less important, since MFV models become more predictive. In MFV, the focus is 
on small, but in many cases visible effects. With increased precision of theoretical methods 
and the amount of data soon available from the LHCb and later, hopefully, from the Super- 
B H , one will have such a level of cross-check among different channels that NP effects of 
even MFV nature should become visible. 

An important test of overall consistency among different flavour processes is certainly 
the global fit to the CKM matrix, which assumes the SM in all observables. However, it 
is also important to consider fits where one restricts to specific classes of observables (in 

1 The same approach, though in a context where effects beyond CMFV are not visible, has been adopted 
in Ref. pi, where the decays Kl — > n°vV, K + — * -k + vv and Kl — > tt°£ + £~ have been studied. 
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particular, suitable ratios thereof) that, when assuming a given NP framework besides the 
SM, are affected in a controlled way or are not affected at all. A first example of this is 
precisely CMFV [|, |, §, that has been analyzed in detail by the UTfit collaboration pQ] . 

However, as first pointed out in Q, the UT analysis suggested by Q does not account 
for the most general framework of MFV. Ref. 0] showed in fact that the definition || of 
MFV does not preclude the ratio AM d /AM s - included in the CMFV UT analyses - to 
be different from the SM value, and therefore, that it should not be included in MFV UT 
fits. 

As a first aim of this paper, we carry out a MFV fit to the UT. When assuming MFV, 
the only quantities allowed to enter the fit are angle determinations and measurements 
of tree-level processes. Concerning the former, we restrict to sin 2(3 as measured from 
B — > ipKs, while we do not include a and 7, for reasons to be explained below. As for 
tree- level processes, we instead restrict to the semileptonic B decays allowing to access \V u b\ 
and \ V c b\. In particular, while the value of \ V c b\ is quite well established, the same cannot yet 
be said about \V u b\, whose inclusive and exclusive determinations are in some disagreement 
with each other, possibly signaling the presence of an underestimated systematic error in 
either of the two determinations. 



It was noted in ref. [11| that, if one uses the inclusive value for \V u b\ in a global UT 
analysis and parameterizes the presence of NP in a model independent way the fit 
adjusts the tension between sin 2(3 and the 'too high' value for \V u b\ by introducing a small 
negative new phase cpB d = —2.9° ± 2.0° in the S^-mixing (see also || |l3|]). 

In the present work, we adopt a somehow complementary point of view, since we focus 
on MFV. The definition of MFV Q precludes the existence of new CP violating phases 
beyond the CKM one. In this framework, we discuss the impact of both the exclusive and 
the inclusive averages for 11/^1 on the UT determination. We show in particular how the 
inclusive | V u b\ favors two non-SM like solutions in the ~p — rj plane. 

Pursuing the possibility that the inclusive \V u b\ determination be the correct one, the 
next question is whether one can find a MFV extension of the SM able to produce either 
of the two non-SM solutions mentioned above. In particular, one can consider the shifts 
in the value of the side Rt implied by the two solutions. Since Rt is (in the SM) related to 
the ratio AM^/ AM S between the B^- and the -B s -system oscillation frequencies, one can 
translate the above shifts into required values for the NP contributions to the ratio itself. 
We consider the explicit example of the MFV MSSM. Our findings indicate that the model 
is able to produce the required amount of corrections only in certain fine-tuned regions of 
the parameter space. This conclusion holds, barring a substantial shift (above 2 standard 
deviations) in the present central values for the low-energy parameters £ and Bk- 



2. MFV fit of the UT 



We now turn to our first task, i.e. the determination of the UT when MFV is assumed. 
In particular we will focus on the UT side Rt, which will be the relevant quantity for our 
subsequent analysis. We have today a whole host of observables, which bear dependence 
on certain combinations of the CKM matrix entries. Hence the determination of the UT 
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rarameter 


Value 


rarameter 


Value 


m c (m c ) 


(1.24 ± 0.037 ± 0.095) GeV 




1.16639 x 10 GeV 


m t [m t ) 


(162.3 ± z.zj Lrev 


JK 


(loy.o ± l.oj MeV 


m K + 


(493.677 ± 0.016) MeV 


B K 


0.79 ±0.04 ±0.09 


Arrix 


(3.4833 ± 0.0066) x 10~ 12 MeV 


a s (m 2 z ) 


0.1176 ±0.0020 unif 


m Bd 


(5.2794 ± 0.0005) GeV 


fs d 


(0.191 ±0.027) GeV 




(5.3696 ± 0.0024) GeV 


B d 


1.37 ±0.14 


mw 


(80.423 ± 0.039) GeV 




1.23 ±0.06 



Table 1: Input parameters used in our fits. Values are taken from the CKMGtter input table 
updated to summer 2006 JT^ /, with the exception of £, for which we take the average performed in 
|T3/. Unless otherwise stated, the Erst error is Gaussian, the second (where applicable) is uniform. 



apex (j),rj) - and all related quantities - follows from a global fit [14, 10 1, which can include 
all or a subset of such observables. 

We consider three fits of the UT, namely 

1. a SM fit, including all the 'classical' constraints, 

2. a MFV fit I, including only tree-level observables and sin2/3^,x S ; 

3. a MFV fit II, analogous to the MFV fit I, but keeping only the inclusive averages for 
both | Kb | and \V cb \. 

Concerning the full SM fit, it can be performed by combining all the available experi- 



mental information (see [14, [lCj]) and assuming the SM. The state-of-the-art results for the 



SM determination of Rt from the CKMfitter and UTfit collaborations (95% CL) read 

CKMfitter: (R t ) S M = 0.868±g;^| , UTfit: (R t )sM = 0.906 ± 0.062 , (2.1) 

in very good agreement with each other. We have performed our SM global fit, using 
the CKMfitter package [14}, p|, a publicly available FORTRAN framework allowing CKM 



analyses in various statistical approaches, e.g. the frequentist one, used in the present 
work. We took advantage of the CKMfitter package also for the other fits performed in 
the present work. It would be interesting to perform an analogous analysis using the UTfit 
code, which is however not yet publicly available. The input parameters and constraints 
used in our global fit is listed in Tables Q and ^. For R t we find (95% CL) 

{R t )sm = 0-8921°;^ , (2.2) 



in good agreement with the findings in eq. (|2.l| ). Eq. (|2.2| ) will be taken as reference figure 
in our subsequent analysis. 

Turning to the MFV fits, the constraints used are collected in Table § and the corre- 
sponding results are displayed in Figs. |l}}2} Some comments are in order on the choices 
of \V u i\. Concerning the inclusive value, we mention that we tried alternatively all the 
averages reported in Ref. [17]. The latter are obtained by analyzing the inclusive data 
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Constraint 



Value 



Ref. SM fit MFV fit I MFV fit II 



\V ud \ 



IK 



us 

a 



7 
AM d 
AM S 

\ € k\ 
sin 2(3 

\Vcb\^ cl) 
\V cb \^ 

|Kfe| (incl) 
|Kf>| (excl) 



0.9738 ± 0.0003 
0.2257 ±0.0021 
(94 ± 8)° 
(83 ± 19)° 
(0.507 ±0.004)/ps 
(17.77 ± 0.12) /ps 
(2.221 ±0.008) x 10~ 3 
0.675 ± 0.026 

(41.7 ±0.7) • 10~ 3 

(39.2 ±0.7 ± 1.4) • 10" 3 

(45.2 ± 1.9 ±2.7) • 10~ 4 

(35.5 ±2.5 ±5.0) • 10" 4 



m 
m 
m 

Hi 





m 




Table 2: Constraints used in the various fits. The first error is Gaussian, the second (where 
applicable) is uniform. A (V ) indicates that the constraint is included in the corresponding ht. 



through use of three alternative theory prescriptions, namely BLNP |]T^, 19, 2C, [H], p2| 
chosen in our final fits, and 



\V ub \^ = (44.6 ± 2.0 ± 2.0) • 10" 4 , HFAG average using DGE || , 
\V ub \( inc V = (50.2 ± 2.6 ± 3.7) • 10" 4 , HFAG average using BLL pi . 



(2.3) 



The choice of the BLNP average is, for our purposes, the most conservative. In fact, in the 
MFV fit II, it leads to two Rt solutions which are closer to the SM one (eq. ( |2.2j )) than in 
the cases where one uses either of the two other averages in eq. Q2.3| ), 

Concerning instead | V^t, | ( cxcl ) 5 other determinations are provided e.g. by Refs. |25] and 
1 26], also quoted in [17|. All of them are consistent with each other and we took the result 
of [^7| for definiteness. 

It can be noted that in the MFV fits of Table || we did not include the constraints on a 



and 7. We mention that inclusion of 7 from the tree-level determination [10] in the MFV 
fit II has the effect of 'selecting' the large- Rt solution (see below) among the two solutions 
displayed in Figs, [l] and g while inclusion of both a and 7 results in a single, SM-like 
solution, with a quite large error. We decided to exclude the a and 7 constraints from our 
main analysis, since our aim here is to make the tension between a 'large' | V^t, | value and 
the present sin 2(3 determination in the context of MFV most transparent. 

The Rt values resulting from the fits of Figs. 1^-0 are the following (95% CL) 



(i^Wv-i G [0.62,0.78] 
ORfWv-lI G [0.54,0.79] 



(#! 2) )mfv-i G [1.07,1.26] , 
(^fWv-ii G [1.06,1.34] . 



These values should be contrasted with the SM result for R t , given in eq. (f 
comments are in order here. 



(2.4) 
(2.5) 

Some 
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Figure 1: Results of the MFV fit I. See Table [| for the constraints used. Left panel displays the 



selected area in the (p — 77) plane, and the UT from the SM fit (2.2 ) is shown as reference. All 
contours are at 95% CL. Right panel reports the corresponding confidence level prohle for R t . 




Figure 2: Results of the MFV fit II. See Table [| for the constraints used. Left and right panels as 
in Fig. 0. 



• In the MFV fit II, the inclusion of only | | C incl ) determines two disjoint solutions 
for R t (eq. ( |2.5| )), which display a discrepancy with respect to the SM solution (eq. 
fl2.2|) ), at more than the 2<r level. 

• The presence of two disjoint solutions for Rt is featured also by the MFV fit I, where 
both the |^ 1 b| ( ' mc ^ and 1 1 C excl ) determinations are included. 

• The two Rt solutions turn out to be even better resolved in the MFV fit I. This is 
due to the inclusion in the fit of the exclusive determinations for both | V u i, | and | V c b \ ■ 
In fact, on the one hand, the latter would shift the central value for the average of 
|Kt&|/|K:b| respectively down and up, as compared to the value obtained in the MFV 
fit II, and the two effects tend to cancel. On the other hand, the presence of two 
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determinations for both \V u b\ and \V c b\ reduces the final [Ktbl/lKbl error with respect 
to the MFV fit II case. 

One possible approach to the above findings is to raise doubts on the | 1 C incl ) deter- 
mination p0| , 31] (see [32] for the exclusive approach and [33] for a very recent reanalysis). 



We observe, however, that the inclusive value for \V u b\ results from an average among quite 
a large number of modes, and the latter display consistency among each other. So the 
average would definitely appear under control, were it not for the above mentioned (and 
well known) discrepancy with the \V u b\ value preferred by the global SM fit, close to the 
exclusive determination. 

In the following we will pursue the possibility that the | V^t, | ( incl ) determination be the 
correct one, and answer the question whether the MFV MSSM can account for either of the 

(1 2) (1 2) 

discrepancies (Rt)sM — (Rt ' )mfv- I n particular, for (R\ ' )mfv we will take the results 
from the MFV fit II, which are the most conservative. Our study addresses the eventuality 
that \V u b\ should consolidate to a value above that required by the SM fit. Since such 
value should most probably lie somewhere in between the present exclusive and inclusive 
determinations, our study represents somehow the 'limiting' case of a 'high' \V u b\ value. 

We finally note that, since the MFV determination of the UT can rely only on a 
handful of reasonably known observables (at present those in Table |2|), our study stresses 
the importance of an accurate determination of \V u b\ and of the angles a and 7. Angle 
determinations are fortunately a major task of the LHCb program and later, hopefully, of 
the Super-B M, which will also allow a precision determination of \V u b\- 



3. Corrections to R f within the MFV MSSM 



3.1 General formula for R t 

In the previous Section, we have focused our attention on the determination of the UT side 
R t . This is because Rt can in turn be expressed in terms of the and B s mass differences 
AMd and AM S . For large tan/3, these quantities can undergo large SUSY corrections even 
in MFV [34]. To display the 'dependence' 2 of Rt on NP contributions to AMd >s , we write 
the following chain of equalities || , starting from the very Rt definition 



Rt 



\VtdVib\ 
\v cd vx\ 



V, 



td 



v, 



cb 



£ frnFs 



AV m B , 



AM, 



(3.1) 



SM 



In presence of NP contributions to AMd )S , one can always write relations of the kind 



AM q = AMg M (l + f q ) , q = d,s 



(3.2) 



where the l.h.s. represents the full theoretical predictions for the mass differences, which 
must in turn be identified with the experimentally measured quantities, since we are sup- 



2 The Rt definition is of course only dependent on CKM matrix entries. 
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posing the presence of NP. Then, eq. (|3.1|) can be rewritten as [||, 



Rt 




where in the last equality we have plugged in the experimental values for AM^g. 

In the light of eq. (3.3), we consider now the two solutions R^ (see eq. (|2.5|)) from 



the MFV fit II of the previous Section. We want to address the question whether the 
discrepancy with respect to the SM solution can be accounted for within the MFV MSSM. 
This will be the case if shifts f q from eq. ( |3.2| ) are sufficiently large to correct eq. ([T^) by 
the required amount. 

3.2 MFV MSSM with low tan/3 

This case, corresponding to tan/? 10, was addressed in detail in Ref. j?]]. There it was 
found that corrections f q are positive and do not exceed a few percent. It was also stressed 
that corrections tend to display alignment between the AM^ and AM S cases, since NP 
contributions from scalar operators able to distinguish between the q = d,s channels come 
with a factor of m^mq/My^ and are below the 1% level. Therefore, within the MFV MSSM 
at low tan/3, the ratio f q = AM^/AM^ tends to be the same for the d, s cases and 
deviations on eq. (|3.3|) are not sufficient to reproduce the Rt solutions of eq. (|2,5| ) . 

3.3 MFV MSSM with large tan /3 

Within the MSSM at large tan /3, the formulae for AM^ and AM S feature also contributions 
from the so-called Higgs double penguins (DP), with exchange of ho,Ho and Aq bosons. 
For values of tan (3 ^ 30, and assuming MFV, the Hq and Aq double penguins usually 
provide the dominant NP contribution which is found to be negative in the whole part of 
the SUSY parameter space with M A G [200, 500] GeV |§, H> 11> 11- The proportionality 



of such contribution to the external quark masses makes it generically negligible for the 
AMrf case, while for AM S sizable corrections are possible. 

In order to address the question whether the SUSY corrections to AM^ and AM S , 
computed in the MFV MSSM at large tan/3, can produce the shifts required for Rt to 



reach any of the solutions in eq. (2J3), we adopt the following strategy 

1. compute the SUSY contributions to AAf^ s in the MSSM with large tan/?; 

2. perform the MFV limit, according to the EFT definition (3|; 

3. study the quantity (1 + /^) /(l + f s ) in the SUSY parameter space left after the MFV 
limit. 

Concerning point 1, the calculation can be accomplished following the procedure of Ref. 
1 34], which allows for a resummation of large tan/3 corrections. For a detailed description 
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of this procedure, we refer to Section 2 of [34|. We mention that in our numerical analysis 



we consistently take into account effects coming from flavour off-diagonal squark mass 
matrices. We do not include, instead, effects arising from large tan j3 corrections to the 
Higgs propagator entering DP contributions. The latter corrections were found to have a 
non-negligible effect only for Ma % 160 GeV |39|. Similar conclusions will be drawn in 
Ref. |40p . Since we do not consider such light pseudoscalar Higgs masses in our numerical 
analysis, these new effects do not change the basic findings of the present work. 

Turning to point 2, the MFV limit was performed in exact analogy with Ref. J?]. In 
particular we expand the soft terms, which enter the squark mass matrices, as functions of 
the SM Yukawa couplings. After such expansions, the soft terms are still parametrically 
dependent on an overall normalization mass scale and on the dimensionless coefficients 
tuning the proportionality to the Yukawa couplings themselves. 

Finally, let us go to point 3. From the above discussion, it is clear that the SUSY 
parameters one needs to consider in the analysis are the same as in the low tan/3 case, 
with the addition of the neutral physical Higgs masses, entering the DP contributions. 
In order to have accurate numerical predictions for such masses, we used the package 
FeynHiggs 41, 42, 44], which calculates the whole physical Higgs spectrum. The only 



additional parameter required by FeynHiggs besides those already present in the low tan (3 
case is Ma- Consequently, in the notation of parameters are 

m , A, Mi , M 2 , M~ g , (3.4) 
H , M A , tan/3 , (3.5) 

with the addition of the 12 'MFV coefficients', governing the proportionality of the soft 
terms to the SM Yukawa couplings (see |]7j). We observe that all these parameters are 
real, since the presence of complex phases would induce non-CKM CP violation in meson 
mixings. The latter would then contradict the MFV hypothesis. 

To explore the SUSY parameter space, we adopt the following strategy. We study the 
ratio (1 + / s )/(l + fd), which enters Rt in eq. (|3.3j), by generating the MFV coefficients 
with flat distributions in the same ranges as (these ranges are also collected in Table |||) . 

Concerning the parameters in eqs. ( |3.4| )-( |3"1)| ), we make the following observation. For 
large tan/3, box contributions are in general much smaller than the DP ones; in addition, 
the 'alignment' of the NP box contributions between the AM^ jS cases, stressed in Section 



3.2, is found to hold for large tan/3 as well. As a consequence, when DP contributions 



are small, one has (1 + /] 3 ° x )/(l + /^ ox ) ~ 1; conversely, when DP contributions are large, 
box corrections are just a small correction, with negligible modifications of the distribution 
shape. 



For the reason above, and as already stressed in Section 3.2, box contributions alone are 
not able to account for the Rt solutions given in eq. (^^) • In order to address the analogous 
possibility for the double penguins, we set m = 1 TeV, thus decoupling box contributions, 
which would partially cancel the (dominant) DP corrections. Parameters A and |^| are 
instead generated flatly in the ranges A £ [—2, +2] TeV and 6 [0.2, 2] TeV. On the sign 
of ji we will come back in the following discussion. The choice of the EW gaugino mass 
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Uniformly distributed 



Fixed 



ai,2,3 G [0.25, 1] 


tan p = {30, 50} 


04,5> fe l,...,8 e [-1, 1] 


M A = {200,500,800} GeV 


A e [-2, 2] TeV 


m = 1 TeV 


\fx\ £ [0.2,2] TeV 


Mi = M 2 = 500 GeV 




M 5 = {200, 500, 800} GeV 



Table 3: Parameter choices made in our scans. For the definition of the MFV coefficients a;, bj, 
as well as of the squark biliear and trilinear mass scales m and A, see Ref. (^j. 



fi>0, tan p = 50 , M A = 200 GeV 



H < , tan p = 30 , M A = 500 GeV 



1200 
1000 
800 
600 
400 
200 



J 



2000 



1500 



1000 



500 



0.4 



0.6 



0.8 



1.2 



0.4 



0.6 



0.8 



1.2 



R, 



Figure 3: Distributions of values for R t , resulting after scanning A, \fi\ and the MFV parameters 
(see text for details). The left panel corresponds to /i > 0, Ma — 200 GeV and tan/3 = 50; the 
right panel corresponds to \l < 0, Ma = 500 GeV and tan/3 = 30. The grey (blue) distributions 
are without (with) the constraint (3.6). The vertical green bands correspond to the Rt solutions 
given in eq. (2.5). 



parameters plays a negligible role and we set them as M\ = Mi = 500 GeV. Finally, we 
considered three discrete values for the gluino mass, namely Mg = {200, 500, 800} GeV. 
Our results turn out to bear negligible dependence also on the choice of the gluino mass: 
we will further comment on this in due course. In our plots, the value M g = 500 GeV is 
chosen for definiteness. 

On the other hand, the quantity (1 + / s )/(l + fd) bears obviously a strong dependence 
on the choice of the mass Ma (which also sets the scale for the H$ DP) and on the 
value of tan /?. We have explored, in discrete steps, the ranges Ma £ [200, 800] GeV and 
tan/3 G {30,50}. 

This completes the description of the choices made for the various SUSY parameters 
in our MonteCarlo study. 3 These choices are also summarized in Table [3|. In Fig. ^ we 

3 We mention that, strictly speaking, our choice of SUSY masses does not always satisfy the condition 
A^susy 2> Ma, required by the effective Lagrangian framework of Refs. p7| . We assume, as also done 
in the existing literature, that the latter approach be applicable also for heavy Higgs masses comparable 
in size with Msusy- A fully rigorous treatment would require a two loop calculation in SUSY, beyond the 
scope of the present work. 
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(i>0, tan (5 = 50, M = 200 GeV 



H<0, tanp = 30, M = 500 GeV 



T 

g"-7 

cc 
m 



-10 




0.5 



1 

AM,/AMf 



CDF + D0 excl. 
(95% CL) 



1.5 




AMj/AMj 



Figure 4: Correlation between BR(B S 
Parameters as in Fig. ||. 



jtx+At") and AM diS in the MFV MSSM at large tan/3. 



show as grey distributions the implied effects on Rt for choices of Ma and tan (3 where such 
effects are significant. The left panel corresponds to the choice Ma = 200 GeV, tan (3 = 50 
and [i > 0, while the right panel corresponds to Ma = 500 GeV, tan/3 = 30 and /x < 0. 

Here we note that, from the left panel of Fig. ||, the long tail looks able to reproduce 
the lower Rt solution given in eq. (|2.5| ). However, a large portion of the tail is actually 
unphysical, since it corresponds to values of AM S which are outside the range allowed by 
the theoretical error (AM^ is only weakly affected by NP contributions, as we saw at the 
beginning of this section). When calculating R t out of AM^ SSM and AM S MSSM , one should 
in fact take into account the constraints coming from the AM q measurements. However, 
these constraints are by far dominated by the corresponding theoretical error. We stress 
that, in our case, the latter must take into account not only the uncertainty coming from 
the lattice parameters, but also the error associated with the CKM entries. For these 
entries one must in fact consistently use the determinations coming from the MFV fit, 
whose associated error is larger than that resulting from the usual SM fit determination. 

The largeness of the final theoretical error makes then the AM q constraints not very 
effective. We decided not to include this 'filter', since a much more effective one is provided 
by the BR(1? S — > fi + [i~) upper bound . The latest combined bound from the CDF and 
D0 collaborations reads Ef 



BR(5 S -> /i + /0 < 5.8 x 10" 



CL] . 



(3.6) 



4 The BR(B S — > [i + fi~) was calculated following exactly the same procedure as the one adopted for 
meson mixings For a recent discussion of these observables using [Q, although in a different context, 
see 



Ref. g|. 

3 For the previous bound from the CDF collaboration, see E^]. 
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The distribution of values for Rt which survives the B s — * fi + jjT constraint is shown in blue 
in Fig. ||. For positive //, the constraint completely excludes the possibility of significant 
corrections to Rt, and sets the latter back to a SM-like value. For negative \x a similar 
conclusion holds, considering that less than 1% of the MC-generated values for Rt lie 
within the 2a band of the lower Rt solution, after the constraint ( |3.6| ) has been taken into 
account. We note that our results in Fig. || can be directly compared with Fig. 3 of p9[ . 
which reports the quantity R sd = (AM s /AMf M )/(AM d /AM^ M ). The latter is related to 
Rt through Rt = 0.913 x \JR sd - One can see that our MonteCarlo allows slightly larger 



maximal effects than in the case of 3J]. This is most likely due to the fact that we do not 



include the constraints from B — > X s j and B — > tv, in constrast to Ref. [39]. Including 
these constraints would however only strenghten our conclusions. 

The severeness of the constraint ( |3.6D on the allowed MFV MSSM corrections to AM S 
is shown in Fig. |j. The latter reports the correlation between BR(i3 s — > and 
AM S in the MFV MSSM 6 at large tan 0, that is roughly given by 

AM S DP oc - (J^J 2 BR(B S - . (3.7) 

An interesting aspect of the above correlation is that, in the regime of DP dominance, 
a large part of the dependence on SUSY parameters other than tan (3 and Ma cancels, 
since it is common to BR(1? S — > n + and AM S . This is in particular true of the 'e- 
factors' ey,e3 and eo (compare eqs. (6.25) and (6.40) of Ref. @]), which incorporate 
the dependence on the SUSY particles entering the loops. This explains, for example, the 
negligible dependence of the correlation ( |3.7D on the choice of the gluino mass, already 
mentioned in the discussion of the parameter ranges for the MonteCarlo. A similar finding 
holds for the ^-distributions of Fig. ^ once the B s — > fi + fJ.~ bound is taken into account. 

In Fig. |], we show two representative cases for the correlation fl3.7|) . As stated pre- 
viously in the text, parameters are chosen in order to have large (correlated) effects. For 
positive [i, Fig. ||| (left panel) shows that the bound ( ^^ ) basically excludes effects on AM S 
exceeding —10%. In addition, the correlation with the prediction on B s —* /i + is largely 
lost in the allowed region. 

For negative fi, neutral Higgs contributions to both BR(£> S —* /i + and AM S are 
generically larger than in the corresponding positive fi case. Furthermore, according to 



eq. (3/7), one can choose larger values for Ma and/or lower values for tan/3 in order to 
have a large correction to AM S while fulfilling the constraint on BR(£? S — > /i + /i~). 7 In 
Fig. |] (right panel) we report a case with these features, corresponding to tan = 30 and 
Ma = 500 GeV. In fact, for negative values of one needs Ma ^ 500 GeV if tan/3 > 30 
in order to observe a significant correlated effect between AM S and BR(i? s — > /i + /i~). One 
should also keep in mind that for fi < the MSSM worsens the (g — 2)^ discrepancy with 
respect to the SM §|, IEI- 

We note here some peculiar features of the MFV MSSM corrections to AM ds , which 
make them unable to reproduce any of the two Rt solutions (^). 



For a discussion of the correlation in the general MSSM we refer the reader to 
7 We thank P. Paradisi for drawing this point to our attention. 
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Figure 5: CL profiles for AM| M (left panel) and AM S SM (right panel) for the MFV fit II. 

a) Among the CKM elements entering the SM formulae for AMd s , only Vtd is signifi- 
cantly modified within the MFV fit with respect to V^ d M . Since Vtd enters only AMd, 
in order to reproduce the experimental values one would need a correspondingly large 
correction fd on AM| M and a neg ligible f 3 correction to AM S SM . The MFV MSSM 
tends instead to give large corrections only to AM S . In fact DP contributions are 
sensitive to the external quark masses, and NP effects scale as m q mb for AM q . 8 

b) In absence of the BR(£> S — ► constraint, corrections to Rt tend to reproduce 
the lower solution, since f s in eq. (|3.3| ) is negative. On the other hand, the higher 
solution looks the favored one, e.g. by the tree-level determination for the angle 7. 

Concerning point a, we show in Fig. | the CL profiles for AMff within the MFV fit II. 
As one can see, in this case the SM formula for AMd favors two solutions, corresponding 
to the shifts in the values of Vtd- On the other hand, AMf M is perfectly compatible with 
the experimental result. 

We have studied the distribution of values for R t as given in eq. (|3.3| ) for different 
choices of the relevant low-energy parameter £, which brings the largest contribution to the 
overall error. 9 Results are displayed in the left panels of Fig. |6|. In addition, when assuming 
MFV, a distribution of values for Rt can be translated into a corresponding distribution 
for 7, with the only additional uncertainty of sin 2/3, to which however 7 is only weakly 
sensitive. Then, similarly to Rt, one can study the dependence of the 7 determination on 
the value of £. Results are shown in the right panels of Fig. ^. The two solutions for 7 

8 In the light of this, a model able to produce large corrections to AMd while leaving those on AM S 
negligible looks quite hard to construct. An interesting case is however that studied by |3jj, flc| ]. 

9 We mention here that the distributions of values for Rt shown in Fig. ^ are calculated for £ set to its 
central value. Taking into account the £ error does not change the relevant features of the distributions. 
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Figure 6: Rt vs. £ and 7 vs. £ plots for fi > (hig-h panels) and \i < (7ow panels). The blue 
bands represent the range of values allowed within the MFV MSSM. The central superimposed light 
blue belts correspond to the present 2o -range for £ from lattice QCD (see Table [l[). The orange 
lines report the SM results. Finally the two horizontal green bands represent the Rt (left panels ) 
or 7 (right panels) solutions from the MFV fit II, eqs. (2.1 ) and (3.8 ). 



corresponding to the Rt determinations of eq. Q2.5| ) are as follows (95% CL) 

(7 (1) )mfv-ii = [21, 48]° , (7 {2) )mfv-ii = [87, 120]° . (3.8) 

They are displayed in Fig. ^ (right) as horizontal green bands. 

We conclude this Section by noting that further support to the above results is provided 
by studying the quantity ex- The latter presents features exactly analogous to AM^. In 
particular: (i) given its dependence on Vfd, in the MFV fits the SM formula gives two 
solutions; (ii) to reproduce the experimental value, one would need relatively large NP 
contributions; (iii) within the MFV MSSM, such contributions are however negligible in 
the full parameter space, so that agreement would again require the parameter Bk to be 
substantially different (> 2a) from the present lattice determination. 



4. Conclusions 

In the present paper we have pointed out that the large value of \V u b\ from inclusive tree- 
level determinations is not only a challenge for the Standard Model, but also for the MSSM 
with MFV. The patterns and the size of modifications of AM^ and AM S in the MSSM with 



-13- 



MFV do not allow a consistent description of the data in this case, unless very significant 
modifications of non-perturbative parameters determinations relative to the lattice ones 
are made. Concerning in particular the modifications allowed to AM S , we conclude that 
scenarios with 200 GeV < Ma < 500 GeV and tan/3 ~ 50 that predict a significant 
suppression for AM S in correlation with an enhancement for BR(i? s — > fi + fi~) have to be 
fine-tuned in order not to violate the new combined bound on the latter decay mode from 
the CDF and D0 collaborations. Relatively large correlated effects can however still occur 
for negative values of /i and large values for Ma > 500 GeV, increasing with increasing 
tan/3 > 30. 

The resolution of the \ V u b\ problem - that is, the discrepancy between the exclusive and 
inclusive determinations of this quantity - calls in the first place for a better theoretical 
control on the relevant hadronic quantities involved in either determination. If the two 
determinations should then agree on a central value in the ballpark of the present inclusive 
average, insight on the tension existing within the SM and the MFV MSSM will require 
a facility like Super-B ||, where also a very precise determination of the angle 7 from 
tree-level processes is possible. As we have seen, the precise knowledge of this angle could 
definitely tell us whether MSSM with MFV has a chance to be a correct description of 
flavour violating processes. 
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